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Abstract. Under natural spectral stability assumptions motivated by previous investigations of the associ- 
ated spectral stability problem, we determine sharp estimates on the linearized solution operator about 
a multidimensional planar periodic wave of a system of conservation laws with viscosity, yielding linearized 
I/^ n — > U' stability for all p > 2 and dimensions d > 1 and nonlinear fl H" — > L*" n H" stability and 
I/^-asymptotic behavior for p > 2 and d > 3. The behavior can in general be rather complicated, involving 
both convective (i.e., wave-like) and diffusive effects. 

1 Introduction 

Nonclassical viscous conservation laws arising in multiphase fluid and solid mechanics ex- 
hibit a rich variety of traveling wave phenomena, including homoclinic (pulse-type) and 
periodic solutions along with the standard heteroclinic (shock, or front-type) solutions. 
Here, we investigate stability of periodic traveling waves: specifically, sufficient conditions 
for stability of the wave. Our main result is to establish bounds on the solution operator 
for the linearized evolution equations, provided that there holds an appropriate spectral 
condition on the linearized operator about the wave. An immediate consequence is that, 
under mild nondegeneracy assumptions motivated by the low-frequency spectral analysis of 
|0Z3] , strong spectral stability in the sense of Schneider \S11iS^S3\l implies linearized and 
nonlinear n LP r\ asymptotic stability, for all p > 2 and dimensions d > 3. 

The one-dimensional study on spectral stability of spatially periodic traveling waves of 
systems of viscous conservation laws was carried out by Oh & Zumbrun [OZlj in the "quasi- 
Hamiltonian" case that the traveling-wave equation possesses an integral of motion, and 
in the general case by Serre |Sel| . An important contribution of Serre was to point out 
a larger connection between the linearized dispersion relation (the function A(^) relating 
spectra to wave number of the linearized operator about the wave) near zero and the 
homogenized system obtained by slow modulation, or WKB, approximation, from which 
the various stability results of |0Z1| . |Sel| may then be deduced. In |0Z3] . we extended 
this important observation of Serre, relating the linearized dispersion relation near zero to a 
multi-dimensional version of the homogenized system. As an immediate corollary, similarly 
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as in [UZlj . [Selj in the one-dimensional case, this yielded as a necessary condition for 
multi-dimensional stability the hyperbohcity of the multi-dimensional homogenized system. 

As noted in |OZ3j , this relation is also the first step in the study of stability and asymptotic 
behavior. In this article, we use the description of low-frequency spectrum carried out in 
|QZ3j to obtain nonlinear stability and asymptotic behavior in dimensions d > 3 by a mod- 
ification of the Bloch decomposition arguments introduced by Schneider in |S1IS2|S3"] : see 
Theorems [H [21 and[3l Here, the main new difficulty is the fact that spectra A(^) bifurcating 
from the translational zero eigenvalue at ^ = are not smooth at the origin, a standard 
feature of hyperbolic-parabolic systems in multi-dimensions |ZS|Z1] . This is not only a 
technical issue, but reflects quite different behaviors in the present vs. previously consid- 
ered cases. Whereas asymptotic behavior in |S1|S2|S3] was purely diffusive, corresponding to 
a Gaussian kernel, the behavior here is convective-diffusive, corresponding asymptotically 
to a convection-diffusion wave in the sense of [HoZlj . Stability and behavior in dimensions 
one and two remain interesting open problems. 

2 Preliminaries 

Consider a system of conservation laws 

ut + ^fHuh,=Y.[B^\u)u,,)^, (2.1) 

3 h k ' 

u G Z^(open) G M", P G M*", B^^ G M*"^", x G M'^, and a periodic traveling wave solution 

u = u{x ■ V — st), (2-2) 
of period X, satisfying the traveling-wave ordinary differential equation (ODE) 

{J2 v,iykB^\u)u')' = (^ u,P{u))' - su' (2.3) 

3,k 3 

with initial conditions 

n(0) = u{X) =: liQ. 
Integrating (|2.3p . we reduce to a first-order profile equation 

VjVkB^^{u)u' = VjP{u) -su-q (2.4) 

3,k j 

encoding the conservative structure of the equations, where 5 is a constant of motion. 
We here briefly review the generic assumptions made in [UZ3|Ser] . Given 
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(j2.4|) admits a unique local solution u{y; a, s, ly, q) such that u(0; a, s, z^, g) = a. Denote by 
X the period, u) := 1/X the frequency and M and -F-^ the averages over the period: 

when ti is a periodic solution of (j2.4|) . Since these quantities are translation invariant, we 
consider the set P of periodic functions u that are solutions of (j2.4p for some triple (s, i^, g), 
and construct the quotient set V := P/TZ under the relation 

(n 7^ (3/i e M; 7; = - /i)). 

We thus have class functions: 

X = X(u), io = f2{u), s = S{u), V = N{u), q = Q{u), M = M(u), = F^{u), (2.5) 

where ii is the equivalence class of translates of different periodic functions. Note that u is a 
nonconstant periodic solution. Without loss of generality, assume S{u) = and N{u) = ei, 
so that ()2.4p takes the form 

B'\u)u' = f\u) - q 
for q = Q{u). Letting X = X[u) and a = u(0) = uq, the map 

(y, a, s, V, q) ^ u{y; o, s, u,q) - a 

is smooth and well-defined in a neighborhood of {X; a, 0, ei,q), and it vanishes at this special 
point. Here and elsewhere, ej denotes the j'th standard Euclidean basis element. We assume: 

(HO) P, B^^ £C^,k> [d/2] + 1. 

(HI) Re a{J2jk ^^j^kB^^) > ^ > for ah e S'^-^ . 

(H2) The map H : M x x R x S'^-"' x M" ^ R" taking {X; a, s, q) ^ u{X; a, s, u,q)-a 
is a submersion at point {X; a, 0, ei,q). 

3 The Evans function 

Without loss of generality taking S{u) = 0, N{u) = ei, u = u{xi) represents a stationary 
solution. Linearizing (j2.ip about u{-), we obtain 

vt = Lv := Y,iB''v,,U - ' (3.1) 

where coefficients 

B^'' := B^^{u), A^v := Df{u)v - {DB^\u)v)u^, (3.2) 
are now periodic functions of xi. 
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Taking the Fourier transform in the transverse coordinate x = (x2, ■ ■ ■ , x^), we obtain 

+ i{^B''^^kvU-i^A^^,v- B^'^k^jv, (3.3) 

where = (^2, ■ ■ ■ ^^d) is the transverse frequency vector. The Laplace transform in time t 
leads us to study the family of eigenvalue equations 

= (L^- - X)w = {B^^w'y - (A^w)' + i Bi^ijw' + B^^ikw)' 

m Ml 
-i^A'^jW- Y B^^ikijW-Xw, (3.4) 

associated with operators and frequency A G C, where '" denotes d/dxi. Clearly, a 
necessary condition for stability of (j2.ip is that (j3.4p have no solutions w for ^ G M'^"^ 
and Re A > 0. For solutions of ()3.4p correspond to normal modes v{x,t) = e^^e^^'^w{xi) of 

(EH). 



Multidimensional stability concerns primarily the behavior of the perturbation of the top 
eigenvalue A = under small perturbations in ^. To study this behavior, we use Floquet's 
theory and an Evans function that not only depends on A but also on ^1 (which corre- 
sponds to the phase shift) and f |G|0Z3| . To define the Evans function, we choose a basis 
{'w^{xi,^, A), . . . , ti;^"(xi, ^, A)} of the kernel of — A, which is analytic in (^, A) and is real 
when A is real, for details see [OZl|Sel] . Now we can define the Evans function by 



(«;')'(X,e,A)-e^^«i(^'y(0,e,A) 



(3.5) 

l<K2n 



where G M.. We remark that D is analytic everywhere, with associated analytic eigen- 
function for 1 < / < 2n. A point A is in the spectrum of if and only if D{X,£^) = 
with^= (6,0- 

Example 31. In the constant-coefficient case, D{X,^) = ^^^"-^(e^'*^''*'^)"^ — e*^^"^), where 
Hi,l = 1, . . . ,2n, denote the roots of the characteristic equation 

^,2^11 + ^(_^i + i ^ ijii^ . + i Y B^'^ik) 

- (i ^ A^ij + Y B'%Cj + XI)\w = 0, (3.6) 

with w = e^^^w. The zero set of D consists of all A and ^1 such that fj.i{X, ^) = 2^i(mod27rz/X) 
for some /. Setting = i^i in (j3.6p . we obtain the dispersion relation {—B^ — iA^ — XI) = 0, 
where A^ = A^^j and B^ = Zj,k B^%(j. 
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4 WKB expansion and the low-frequency limit 



We now recall the results of [0Z3| describing low-frequency spectral behavior. As a conse- 
quence of (HO), (H2), there is a smooth n + d dimensional manifold V of periodic solutions ii 
in the vicinity of u, where d is the spatial dimension. On this set, one may obtain, rescaling 
by (x, t) — > (ex, et) and carrying out a formal WKB expansion as e — > a closed system of 
n + d averaged, or homogenized, equations 

j 

^ti^2N{u)) + V.,{^2S{u)) = (4.1) 

in the (n + d)-dimensional unknown ii, expected to correspond to large time-space behavior, 
with an additional constraint 

curl (QN) = (4.2) 

coming from the derivation of the formal expansion: specifically, the assumption that QN 
represent the gradient Vx4> of a certain phase function (j){x,t). Here, i7, M, etc. are defined 
as in (j2.5p : see |0Z3| for details. 

The long-time behavior of perturbations of u can thus be studied formally by considering 
the linearized equations of (j4.ip about the constant solution u{x,t) = , u, yielding 
the homogeneous degree n-\- d linearized dispersion relation 

Z^(C,A):=det(A^^— ^(^z) + 5^.e,^^— ^(n)l =0. (4.3) 

Alternatively, it may be studied rigorously through low-frequency expansion of the Evans 
function D(^,A), ^ G M*^, A G C, which yields after some standard but somewhat involved 
manipulations 

Z)(C,A) = Z\i(C,A) + 0(|e,Ar2), (4.4) 

where A\ is a homogeneous degree n + 1 polynomial expressed as the determinant of a 
rather complicated 2n x 2n matrix in (^, A). The zero set (C,A(^)) of A\ thus determines 
the linearized dispersion relation for (j2.ip . with A(^) running over the tangent cone at ^ = 
to the surface of low- frequency spectrum of L as ^ runs over M*^. 

Our main result in |0Z3| was the following proposition relating these two expansions, gen- 
eralizing the result of |Sel] in the one-dimensional case. Define 

Z\(e,A) :=Ai-'^i(e,A), (4.5) 

where A is defined as in (14.31). 



Proposition 4.1. lOZSf Under assumptions (H0)-(H3), Ai = FqA, i.e., 

Z)(C,A) = roZ^(C,A) + 0(le,Ar+2) (4.6) 
Iq 7^ constant, for |^,A| sufficiently small. 
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That is, up to an additional factor of A'^^^ (corresponding to spurious modes not satisfying 
constraint (j4.2t see |QZ3j for further discussion) the dispersion relation (j4.3p for the averaged 
system (|4.ip indeed describes the low-frequency limit of the exact linearized dispersion 
relation D{^,X) = 0. 

Corollary 4.1. JOZSj/ Assuming (H0)-(H3) and the nondegeneracy condition 



det 



(M,.,).o, 



then for X,S^ sufficiently small, the zero-set ofD{-,-), corresponding to spectra ofL, consists 
o/ n + 1 characteristic surfaces: 

A,(0 = -ia,(e)+o(e), i = l,...,n + l, (4.8) 

where aj{(,) denote the homogeneous, degree one eigenvalues of 

excluding {d — 1) identically zero eigenvalues associated with modes not satisfying \4.^ . 

Corollary 4.2. JOZSj/ Assuming (H0)-(H3) and the nondegeneracy condition ( |^. 7| j, a nec- 
essary condition for low-frequency spectral stability ofu, defined as Re X < for D(^, A) = 0, 
^ G W^, and \C, X\ sufficiently small, is that the averaged system be "weakly hyperbolic" 
in the sense that the eigenvalues of A{^) are real for all ^ G M.'^. 

A consequence of Corollary 14.11 is that Xj{(,) are differentiable in |^| at the origin for fixed 
angle but in general have a conical singularity in ^ at = 0, since the eigenvalues 
aj(^) of first-order system (j4.ip are homogeneous degree one but typically not linear. The 
low-frequency expansion of Corollarv 14.11 substitutes in our analysis for the usual spectral 
perturbation analysis by formal series expansion/Fredholm alternative in the standard case 
(as in |S1|S2|S3"] ) that Aj(^) vary smoothly in 

Remark 4.1. The low- frequency stability condition of Corollary 14.21 has been verified 
numerically for the example of isentropic van der Waals gas dynamics |0j , for which periodic 
solutions are known to appear. On the other hand, it was shown in |0Z1| that high-frequency 
instabilities appear for these waves, so they are not in the end stable. 

5 Bloch— Fourier decomposition and the spectral stability conditions 

Following |G|S1|S2|S3| . we define the family of operators 

= e-*«i^iL^-e^«i^i (5.1) 



6 



operating on the class of periodic functions on [0,X]; the (L^) spectrum of is equal 
to the union of the spectra of all with real with associated eigenfunctions 

w{xi,l\) :=e^«i"ig(xi,6,e,A), (5.2) 

where q, periodic, is an eigenfunction of L^. By continuity of spectrum, and discreteness of 
the spectrum of the elliptic operators on the compact domain [0, X], we have that the 
spectra of may be described as the union of countably many continuous surfaces Aj(^). 

The spectrum of each may alternatively be characterized as the zero set for fixed ^ of 
the periodic Evans function Z)(^i, ^, A); see [G|0Z3) . More, a fundamental result of Gardner 
[G] is that the order of vanishing of the Evans function in A is equal to the multiplicity of 
A as an eigenvalue of L^. Thus, we have a description of the eigenstructure of through 
Corollary 14.11 for A, ^ sufficiently small in terms of the characteristics of the first-order 
hyperbolic system (j4.ip . 

Without loss of generality taking X = 1, recall now the Bloch-Fourier representation 

<^) = (irY r [ e'^-^u{i,xi)diidi (5.3) 

of an function u, where u{£^,xi) := ^j^, e^'^*'^^iu(^i + 27rA;,^) are periodic functions of 
period X = 1, n(^) denoting with slight abuse of notation the Fourier transform of u in the 
full variable x. By Parseval's identity, the Bloch-Fourier transform u{x) u{S,,xi) is an 
isometry in L^: 

\Ml^{x) = \\u\\l^(^;L^(xi)), (5-4) 

where -L^(xi) is taken on [0, 1] and L?'{C) on [— 7r,7r] x M'^"^. Moreover, it diagonalizes the 
periodic-coefficient operator L, yielding the inverse Bloch-Fourier transform representation 

e'-^^^o = i^Y r [ e'^-^e'^^'uo{C,x,)dCidi (5.5) 
V27r/ J_^J^d-i 

relating behavior of the linearized system to that of the diagonal operators L^. 

Along with the assumptions (H0)-(H2), we assume the strong spectral stability conditions: 

(Dl) a{L) C {ReA < 0} U {0}. 

(D2) Rea{L^) < -6l[^|2, 6* > 0, for ^ € M.'^ and |^| sufficiently smaU. 

By Corollary 14.1^ either of (Dl), (D2) implies that the eigenvalues aj(^) of (14. 9p are real. 
We make the further nondegeneracy hypothesis: 

(H3) The values aj(^) are distinct. 

Conditions (D1)-(D2) are exactly the spectral assumptions of |S1|S2|S3"] . corresponding to 
"dissipativity" of the large-time behavior of the linearized system. Condition (H3) corre- 
sponds to strict hyperbolicity of the averaged system (j4.ip : presumably it could be removed 
with further effort/alternative hypotheses. 
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6 Linearized stability estimates 



By standard spectral perturbation theory [K], the total eigenprojection -P(^) onto the 
eigenspace of associated with the eigenvalues Aj(^), j = 1, . . . , n + 1 described in Corol- 
lary UTT] is well-defined and analytic in ^ for ^ sufficiently small, since these (by discreteness 
of the spectra of L^) are separated at = from the rest of the spectrum of Lq. Introduc- 
ing a smooth cutoff function 0(^) that is identically one for [^| < e and identically zero for 
> 2e, e > sufficiently small, we split the solution operator S{t) := e^* into low- and 
high-frequency parts 

S'{t)uo := r [ e'^-^mP{Oe''^'M^,xi)d^,di (6.1) 

^^^^ J-TrJR'i-^ 

S"{t)uo := (^Y r [ e^f-^(I-</)P(0)e^«*tio(e,xi)d6<ie (6.2) 



and 



6.1 High-frequency bounds 

By standard sectorial bounds |He|Pa| and spectral separation of Aj(^) from the remaining 
spectra of L^, we have trivially the exponential decay bounds 

||e^«(/-,/<P(0)/||L2([o,x])<C7e''*!l/!lL2([o,x]), (6.3) 

\\e'^Hl-'PPimccf\\LH[o,x])<Ct~h'''\\fh2^[o,x]), (6.4) 

Wd.e'^^il - ^PiO)f\\L^i[o,x]) <Ct-h-''\\f\\m[o,x]), (6.5) 
for 0, C > 0. Together with (|5.4p . these give immediately the following estimates. 

Proposition 6.1. Under assumptions (H0)-(H3) and (D1)~(D2), for some 9, C > 0, and 
all t > 0, 2 < p < CO, 

\\S'\t)f\\L^,)<Ce-''\\fh2(^,), (6.6) 
\\d.S'\t)f\\m,^, \\S''mJh2^,) < Ct-5e-^*||/||i2(,)), (6.7) 

< Ct-^(^-i)e-^*||/||i2(,). (6.8) 

Proof. The first two inequalities follow immediately by (|5.4p . The third follows for p = oo 
by Sobolev embedding from 

and 

which follow by an application of (j5.4p in the xi variable and the Hausdorff- Young inequality 
II/IIl°°(xO ^ variable x. The result for general 2 < p < oo then follows by 

interpolation. ■ 
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6.2 Low- frequency bounds 



Denote by 

G'{x,t;y):=S'{t)6y{x) (6.9) 
the Green kernel associated with , and 

[G|(xi,t;yi)] := </.(e)P(0e^«*[<5s/i (^i)] (6-10) 

the corresponding kernel appearing within the Bloch-Fourier representation of , where 
the brackets on [G^] and [6y] denote the periodic extensions of these functions onto the 
whole line. Then, we have the following descriptions of G^ , [G|], deriving from the Evans 
function analysis of Corollarv 14. 1[ 

Proposition 6.2. Under assumptions (H0)-(H3) and (D1)-(D2), 

n+l 

[G|(xi, t; yi)] = m e^^^^^'qji^, xi)g, (C, yi)* , (6.11) 
i=i 

G'ix,t;y) = ^l-Y J^^e^i-(-y)[Glix,,t;y,M (6-12) 

= (_\ / e*«-(-^),/.(0 5]e^^(5)%(e,xi)g,(^,yi)*d^, 
Vivr/ J^d ^ 

where * denotes matrix adjoint, or complex conjugate transpose, qj{^,-) '^''^d qj{£,,-) are 
right and left eigenf unctions of associated with eigenvalues Aj(^) defined in Corollary 
\4-l\ normalized so that {qj,qj) = 1, where Aj/|.^| is a smooth function of \(,\ and ^ := ^/|^| 
and Qj and qj are smooth functions of \^\, ^ := ^/|^|, and xi or yi, with 5)?Aj(^) < — 

Proof. Smooth dependence of \j and of q, q as functions in L^[0, X] follow from standard 
spectral perturbation theory ^ using the fact that Aj split to first order in |^| as ^ is varied 
along rays through the origin, and that varies smoothly with angle ^. Smoothness of qj, qj 
in xi, yi then follow from the fact that they satisfy the eigenvalue equation for L^, which has 
smooth, periodic coefficients. Likewise, ()6.1ip is immediate from the spectral decomposition 
of elliptic operators on finite domains. Substituting (|6.9|) into (|6.ip and computing 

8y{^^ xi) = e^^'^'^^Syii + 27rA:ei) = ^ ^2^ikx^ ^-tt,-y-2^tky^ ^ e-^^-y[5y{x)], (6.13) 

k k 

where the second and third equalities follow from the fact that the Fourier transform either 
continuous or discrete of the delta-function is unity, we obtain 

G\x,t;y) = i^Y r [ e'^--<pPiC)e'^^'6y{txi)d^ (6.14) 

= (^) VI 

yielding (I6.12P by (|6.10p and the fact that (j) is supported on [— 7r,7r]. ■ 
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Proposition 6.3. Under assumptions (H0)-(H3) and (D1)-(D3), 

sup||G^(-,i,;y)||iP(,), sup ||5,,j,G^(-, t, ; < C(l + (6.15) 

y y 

for all 2 < p < oo, t > 0, where C > is independent of p. 

Proof. From representation (|6.11|) (ii) and 5RAj(^) < — we obtain by the triangle 
inequality 

<C||e-^'«''V(e)llLi(C) <C^(l+i)"^ (6-16) 

verifying the bounds for p = oo. Derivative bounds follow similarly, since derivatives falling 
on Qj or ijj are harmless, whereas derivatives falling on e*^'^^~^^ bring down a factor of ^, 
again harmless because of the cutoff function cp. 

To obtain bounds for p = 2, we note that (|6.11l fii) may be viewed itself as a Bloch-Fourier 
decomposition with respect to variable z := x—y, with y appearing as a parameter. Recalling 
(j5.4p . we may thus estimate 

sup||G^(x,t;?/)||i2(^) = Vsup||0(C)e-^^(^)*gj-(-,2;i)Q'*(-,yi)||i2(^.i2(^^g[o,x])) (6-17) 
y J y 

< C V sup ||0(Oe"^'^''* IIl2(0 ll^j IIl2(o,x) ll^j IIl°°(o,x) (6-18) 
<C{l + t)-i, (6.19) 

where we have used in a crucial way the boundedness of qj] derivative bounds follow simi- 
larly. 

Finally, bounds for 2 < p < oo follow by L^'-interpolation. ■ 

Remark 6.1. In obtaining the key L^-estimate, we have used in an essential way the 
periodic structure of qj, qj. For, viewing as a general pseudo differential expression rather 
than a Bloch-Fourier decomposition, we find that the smoothness of qj, qj is not sufficient 
to apply standard — > bounds of Hormander, which require blowup in ^ derivatives at 
less than the critical rate found here; see, e.g., [H] for further discussion. 

Remark 6.2. Computation (j6.13p applied to the full solution formula (|5.5p yields the 
fundamental relation 

G{x,t;y) = (^Y r I e'^<-'y\G^{xi,t-y^m (6-20) 
\2ttJ J_^J^d-i 

which, provided cr{L^) is semisimple, yields the simple formula 

Gix,t;y) = £ ^^_^e^«-(--^)5^e^^«)*(?,(e,xi)^,(e,yi)*d^ 
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resembling that of the constant-coefficient case, where Xj runs through the spectrum of 
Lg. Relation (j6.20p underlies both the present analysis and the technically rather different 
approach of [0Z2j , with the basic idea in both cases being to separate off the principal part 
of the series involving small Aj(^) and estimate the remainder as a faster-decaying residual. 

Corollary 6.1. Under assumptions (H0)-(H3) and (D1)~(D2), for all p>2, t> 0, 

\\S\t)f\\Lv, \\d,S\t)f\\Lv, \\S\t)dJ\\Lv <C{l + tri^^-l^\\f\\Li. (6.21) 

Proof. Immediate, from ()6.15p and the triangle inequality, as, for example, 
\\S\t)f{-)\\Lv= [ G\x,t;y)f{y)dy _ < I sup t; 



LP{x) 



y 



Additional estimates. For general interest, we state also some easy-to-obtain generaliza- 
tions, even though they are not needed in the analysis. From boundedness of the spectral 
projections Pj{C) = qjiqj,-) in L^[0,X] and their derivatives, another consequence of first- 
order splitting of eigenvalues Xj{C) at the origin, we obtain boundedness of (t){i)P{C)^^^^ 
and thus, by (j5.4p . the global bounds 



\\S\t)f\\L2(,), \\d^S\t)f\\L2(,), < for all t > 0. (6.22) 

By Riesz-Thorin interpolation between (j6.22p and (j6.2ip . we obtain the following, appar- 
ently sharp bounds between various L"? and U*. 

Corollary 6.2. Assuming (H0)-(H3) and (D1)-(D2), for all I < q < 2, p > 2, t > 0, 

\\S\t)f\\LP, \\d,S\t)f\\L,, \\S'it)dJ\\LP<C{l + tr^^-^'-v^\\f\\L,. (6.23) 



6.3 Short-time bounds 

Finally, we recall the following short-time bounds on the full solution operator S{t) = e^*, 
following from standard semigroup theory [Pa] for the second-order elliptic operator L. 

Proposition 6.4. Assuming (H0)-(H3) and (D1)-(D2), for I < p < oo, < t < 1, 

\\S{t)f\\LP<C\\f\\LP, (6.24) 
\\d.Sit)f\\LP, \\Sit)d,f\\LP < Ct~^\f\\LP. (6.25) 
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6.4 Linearized stability in dimensions d > 1 

Theorem 1. Assuming (H0)-(H3), spectral stability (D1)-(D2), implies D 
asymptotic stability of the linear equation \2.1\) . for all p > 2 and dimensions d > 1, with 

\\Sit)uo\\L^' <C{l + ty^^^-p^\\uo\\LinLP, (6.26) 
\\S{t)d,cUo\\LP < Ct"5(l + t)^-|(i-?)||uo||^i^^p for all t > 0. (6.27) 

Proof. Immediate, from ([O]) . (jOT]) . and ([QI]) . ■ 

7 Nonlinear stability in dimensions d > 3 

Define now the perturbation variable v := u — u for u a solution of (j2.ip . 

Proposition 7.1. Assuming (H0)-(H3), letvo eH^,k> [d/2] + l as in (HO), and suppose 
that for < t < T, the norm of v remains bounded by a sufficiently small constant. 
There are then constants 9i^2 > so that, for all <t <T, 

||^;(t)|ll^. <Ce-^^*||^;(0)||^, +C f e-'^^'~'Mh{s)ds. (7.1) 

^0 

Proof. Subtracting the equations for u and u, we may write the nonlinear perturbation 
equation as 

vt + ^{dfj{u)v)^^ - ^{Bjk{u)v^^)xy, = ^{{Bjk{u + v)- Bjkiu))u^^)x,, (7.2) 

j j,k j,k 

- ^ifji^ + V)- fj{u) - dfj{u)v)a:y 
j 

In the uniformly elliptic case 

jk 

we may take the inner product in x of X^^^g ^x'^^ against (|7.2|) . integrate by parts, apply 
Garding's inequality, and rearrange the resulting terms, to arrive at the inequality 

dt\\v\\l.it)<-e\\d^,+\\\l, + c\\v\\l„ 

where 6 > 0, for some sufficiently large C > 0, so long as ||f ||//fe remains bounded. Using 
the Sobolev interpolation ||f ||^fe < C'~^||9^+^f ||^2 + C'||'y||/,2 for (7 > sufficiently large, we 
obtain Ot||u|||^fe(t) < —^11^11^^ +C|[?;|[^25 from which (|7.ip follows by Gronwall's inequality. 
In the general case, we perform an analogous pseudodifferential estimate using a frequency- 
dependent symmetrizer to obtain the same result. We omit the (standard) details. ■ 
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Theorem 2. Assume (H0)-(H3). Then, spectral stability, (D1)-(D3), implies nonlinear 
asymptotic stability of u from L} n H'^ L^, k > [d/2] +1 as in (HO) and p > 2, for 
dimensions d>2, with 

\\u{-,t)-u\\Lv < C{1 + ty^^'^-p'>\\uo - uW^inH^, (7.3) 
\\u{-,t) -u\\jjk <C{l + t)~i\\uo -uW^inHk for alH > 0. (7.4) 

Proof. Taylor expanding about u, we obtain the alternative perturbation equation 

vt-Lv = Y,Q'{v,'^v),^, (7.5) 
j 

where 

Q^v,Vv) = 0{\v\'^ + \v\\Vv\) (7.6) 
as long as \v\ remains bounded by some fixed constant. By Duhamel's principle, we have 

vi;t)=S{t)vo+ f S{t-s)y^dy^Q^{-,s)ds. (7.7) 
Case n H'^. Define now 

r]{t) :=supo<,<J|t;(-,s)||Lp(l + s)l (7.8) 

By Proposition 17. H ||v||j:/2, hence ||i't||2,2, remains small so long as rj remains small, hence 
r] remains continuous so long as it remains small. We first establish 

r]{t) < C{r]o + vit?), m ■= lbo||LinJyfc> (7-9) 

from which it follows by continuous induction that ri{t) < 2CrjQ for t > 0, if 7?o < 1/4C, 
yielding by ()7.8p the result (j7.3p for p = 2. This in turn yields (j7.4p by Proposition 17.11 

By Proposition 17. H (j7.8p . and Sobolev embedding, 

\Q'{v,Vv){-,t)\L^r.Lv < {\Vv\L, + \v\L2)\v\L2nL^ <Cvitf{l + ty-^ (7.10) 

for all 2 < p < oo, in particular p = 2. Substituting into (|7.7p and using (|6.26p . (j6.27p . we 
thus obtain 

\v{;t)\L2<C7jo{l + tri (7.11) 
* ... 



+ C7]{tf / (l+i-s)2-4(t-s)-2(l+s) 

Jo 

<C{r,o + 7i{t?){l + tr^''^-i'^~'f^ 



2 



so long as (1 + s) 2 is integrable, i.e., for d > 3. Noting that 1 — ^ < —j for d > 2, we 
obtain (|7.9p as claimed. This completes the proof of (|7.3p - (|7.4p for p = 2. 
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Case LP, 2 < p < oo. Substituting again into (|7.7p and using (|6.26|) . (|6.27|) . we obtain 

\vi;t)\L2<Crjoil + t)-^^'-'p^ (7.12) 

+ C'n{tf / (i + t-s)5-l(t-s)"^(i-i)(i + s)-| 



Jo 

for d >3. This completes the proof of (|7.3|) for 2 < p < oo. ■ 

Remark 7.1. In the critical dimension d = 2, our stability argument barely fails, due to 
the appearance of a logt term coming from the expression /(I + s)~'^ds. 



8 Asymptotic behavior in dimensions d > 3 
8.1 Second-order approximation 

Expressing ^ in polar coordinates (r, where ^ = r^, r = |^|, ^ = i^/li^l, and expanding 

^3ir,i)=a,{i)r + bjiiy + 0{r'), 

define now 

X*{0-=a,iOr + b,iiy, j = l,...,n + l, (8.1) 

n+l 

[Gt(xi,t;yi)] :=<A(OE^^^'^^^*«j((0,|),xi)g,((0,O,yi)*, (8.2) 

Gt(x,t;y) := (i-)' e*€-(-^)[Gt(xi, t; yi)]de, (8.3) 

noting that q{{r,S,),xi) and i?((r, ^),yi) by Proposition 16.21 are smooth in all arguments. 
Define likewise 

v\x,t):= f G\x,t;y)voiy)dy= [ G\x,t;y){uo - u){y)dy. (8.4) 



Proposition 8.1. Assuming (H0)-(H3) and (D1)-(D3), for all p > 2, d > 3, t > 0, 

\\u{;t)-U-V^{;t)\\LP < C {I + t)'^'^^' \\uo - u\\ L^nH^ , (8-5) 

\\u{-,t) - u- v\-,t)\\Hk < C{l + t)~i~^\\uo - ulliinHfe- (8-6) 

Proof. By smoothness of qj, qj in angle ^, 

\qj{{r,i),xi)-qj{{0,i),xi)\, |^i((r, 0, a;i) - g,((0, 1), xi)| = 0(|C|), 

whence we easily obtain (|8.5p and ()8.6p by estimates like those used in the proof of Theorem 

[2j We omit the details; see, for example, [HoZl| for similar computations. ■ 
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8.2 Reduction to constant-coefficients 

Noting that Span{gj}(0, ^) and Span{q^j}(0, ^) are independent of tlie angle ^, corresponding 
for eacli ^ to the right and left zero eigenspaces and IJq of Lq, we may choose a fixed 
real (since Lq has real coefficients and eigenvalue is real) pair of dual bases tTj and ttj of 
and IJq, (vTj, vTfc) = 6j , j = 1, ... ,n + 1, to obtain 

q,iiO,i),xi) = ^akjiO,i)7Tkixi), (8.7) 

k 
k 

(8.8) 

for some smooth nonsingular matrix- valued functions a, a £ with 

a*a = In+i- 

Denoting by aj, aj the jth columns of a, a, we obtain the Floquet-type factorization 

Gt(x, t- y) = n{xi)g^x - y, t)n{yiY, (8.9) 
g\x-y,t):=(^Y [ e^^<^-y^[gUtm, (8.10) 

n+l 

[glit)] ■■= mT.e^^^^K,{o,i)a,{o,ir, (8.11) 

where 77 := (vri, . . . , iTn+i), fl := (vri, . . . , tTu+i), and 51 is a constant-coefficient operator in 
the sense that it is invariant under spatial translations. Recall, by definition, that 77*77 = 7„, 
so that 77* is something like a pseudoinverse of 77 (not "the" pseudoinverse, however, except 
in the self-adjoint case Sq = IJq). 

Indeed, we may factor a bit further, as 

g'f = W*K = K*W, (8.12) 

where 

Wiz,t):=(—) / e^«-cA(OEe"'^^^*"j(0'^>^(0'^~)*^? (8-13) 
denotes the hyperbolic part of the solution operator and 

7^(z,t):=(-) / e*«-^<A(6E^''^^^*«^(0'^>^(0'^")*^^ (^.14) 
Vivr/ J^d ^ 

denotes the diffusive part, with hj{^ := \^^hj(X). This gives a description of the constant- 
coefficient solution operator s^{t)f := g^{-,t) * f as the composition = SwSk = SrSw 
of commuting hyperbolic and parabolic solution operators Sw{t)f '■= W{-,t)*f, Sxit)!' '■= 
K{-,t)*f, and of the Green kernel g^ as a linear convection-diffusion wave W*K = Sw{t)K 
as defined in |HoZl] for general constant-coefficient systems. 
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8.3 Convergence to linear convection diffusion wave 

Define now w{x,t) := n*{xi)v^ {x,t) G M"+^. Evidently, 

{x,t) = n{xi)w{x,t), w = g'^*WQ, wo{x) := n*{xi)vo{x). (8.15) 

Denote by 

W := [ wo{x)dx= [ n{xi)*voix)dx (8.16) 

the total mass of wq. 

Theorem 3. Assuming (H0)-(H3) and (D1)-(D3), ford>3,t> I, 

\\u{-,t) -u- ngH-,t)W\\Lv < Cil + ty^^^-i^^-^WvohinH" + I^^^oIIlO, (8.17) 
\\u{;t) -u- ng^{-,t)W\\H, < C{l + tyi''2{\\vo\\LinH>' + \\xwo\\li), (8.18) 

where W G M""'"^ is the constant mass vector defined in i8.16]) . Moreover, 

C,il + t)'i < \\gH;t)\\L2<C2il + t)-i, (8.19) 

so that in general ||n(-, t) — u — ng\-,t)W\\i2 « WlJg'^ {■,t)W\\i2 , or u{-, t) — u IJg^W. 

Proof. By Proposition 18. 11 (I8.15p . and boundedness of il and derivatives, it is sufficient to 
show that 

y * wo - g^W\\LP < Cil + t)-5(i-i)-5 WxwoWli, 
\\g^ *Wo - g^WWijk < C(l +t)"^"5||a;u;o||Li, 

or, by Hausdorff- Young's inequality and Parseval's identity, that 

y{wo{0-W)\\L.<C{l+t)--^^^-"p^-"^\\xwo\\Li, (8.20) 
11(1 + |ehP(u'o(e) - VF)||i2(^) < C{l+t)-l-^\\xwo\\L^, (8.21) 



that 



where ^ + ^ = 1- Noting, by the Mean Value Theorem and Hausdorff- Young's inequality. 



||wo(C)-W^IIl°=(0 = \\wo{0-wo{0)\\l^{^) < \^\\\d^wo\\L^{^) < |^|||0^wo||l°o(o < ICIIkw^ollLi 
we readily obtain (f8:20]) - ([8:2T]) from 



\m\L.<c{i+t)-^2('--p'^-'2 and iie(i + ieiViiL2(.)<c(i+o-^^ 



as follow by direct computation from representation (|8.10p - (18.1ip as, likewise, does (j8.19p : 
see [HoZl] for similar computations. ■ 
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9 Discussion and open problems 



Theorem [3] shows that the L-^-asymptotic behavior of perturbations v with initial value vq 
possessing an first moment is given by a linear convection-diffusion wave 



with amplitude determined by the modulated mass M = J^a n* {xi)vo{x)dx, decaying in 
at the rate of a heat kernel (Gaussian). Due to the convective-diffusive structure of 
, this is essentially all we can say without further assumptions on the explicit structure 
of the hyperbolic system (j4.ip . For, geometric effects such as focusing or defocusing of 
characteristics can greatly affect the norm W * K for norms p > 2, as discussed in 
[HoZlj for the specific case of the wave equation. It might even be that for sufficiently large 
p a different part of the solution dominates behavior. 

A brief consideration reveals that the zero eigenspace of Lq spaned by the columns of IJ 
consists of tangent directions along the manifold of possible periodic solutions nearby u. 
Thus, our description (|8.9|) of lowest-order behavior as the product of 77 and a solution 
g"^ * {ttvq) of a diffusive regularization of a hyperbolic system corresponding to (j4.ip can be 
viewed roughly as a linearized version of the formal description by WKB approximation, in 
which, to lowest order, the solution is approximated by 



where C indexes the manifold of nearby traveling waves and ip is a scalar phase function 
with V^V = Nf2; see |Sel|UZ3] for further discussion. 

For lower dimensions d = 1 and 2 where decay of the linearized solution is slower, behavior 
is not expected to be dominated by its linear part, and indeed we have seen that the 
description of the solution as linear part plus error is too crude even to close a stability 
argument. A very interesting direction for further investigation would be to attempt to 
encode the lowest-order behavior at a nonlinear level using (j9.ip or a slight modification, so 
as to eliminate the largest terms in the nonlinear residual and close the stability argument. 
See, for example, the argument used in jHoZ2| to obtain stability and behavior of scalar 
viscous shock fronts in the critical dimension d = 2, in which quite similar difficulties arise. 
See also the remarkable work of Schneider |S1|S2|S3] on stability of patterns in dimensions 
d = 1,2, in which the nonlinear behavior encoded by modulation equations likewise plays a 
crucial role in the analysis by revealing unexpected cancellation needed to close the stability 
estimates. 

We remark that our way of getting linearized estimates based on isometry properties is 
essentially different from, and somewhat simpler than, either the weighted norm approach 
of |S1|S2|S3] or the one-dimensional pointwise approach of [QZlj . This more primitive ap- 
proach allows us to treat cases, as here, for which the low-frequency dispersion relation is 
not smooth, as can be expected for general systems for which convection plays a role. 



n{xi)g\x, t)M = n{xi)W * K{x, t)M 



u 
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Finally, we recall that the analysis of [S1IS2|S3] concerns general multiply periodic waves, 
i.e., waves that are either periodic or else constant in each coordinate direction. It would 
be very interesting to consider whether such waves with two or more periods can arise as 
solutions of conservation laws, and if so, what would be the resulting behavior, even at a 
formal WKB level. 
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